Relativistic stars in f(R) gravity 
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We study the strong gravity regime in viable models of so-called f[R) gravity that account for 
the observed cosmic acceleration. In contrast with recent works suggesting that very relativistic 
stars might not exist in these models, we find numerical solutions corresponding to static star 
configurations with a strong gravitational field. The choice of the equation of state for the star 
is crucial for the existence of solutions. Indeed, if the pressure exceeds one third of the energy 
density in a large part of the star, static configurations do not exist. In our analysis, we use a 
polytropic equation of state, which is not plagued with this problem and, moreover, provides a 
better approximation for a realistic neutron star. 
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One of the most challenging tasks for cosmology and 
fundamental physics today is to try to understand the 
apparent acceleration of the Universe. Among the many 
models which have been proposed, a class which has at- 
tracted a lot of attention is the so-called f{R) gravity the- 
ories where the standard Einstein-Hilbert gravitational 
Lagrangian, proportional to the scalar curvature R, is 
replaced by a function of R while the matter part of the 
Lagrangian is left unchanged (see e.g. [l| for a recent 
review). After many detours, it has been realized that 
viable f{R) theories must satisfy stringent conditions in 
order to avoid instabilities and to be compatible with 
the present laboratory and astrophysical constraints. A 
few models (H-Q have been carefully constructed to meet 
these requirements, using in particular the chameleon 
mechanism [5| to satisfy the solar system constraints (the 
binary pulsar constraints are, so far, weaker for these 
models Q). 

However, recent works 0, Q studying the strong grav- 
ity limit of these models have questioned their viability 
by suggesting that, inside neutron stars, where the effects 
of general relativity are strong, one would easily reach 
the singularity where R becomes infinite. In particular, 
in this impression seemed reinforced by the impossi- 
bility to construct numerically relativistic stars beyond 
some critical value of the gravitational field. 

In this Letter, we reexamine this question and, unlike 
these previous works, find that highly relativistic stars 
can be obtained numerically. We also give qualitative ar- 
guments to understand the existence of these solutions. 
An important conclusion of our work is that the equation 
of state must satisfy p—3P > in most of the star. Oth- 
erwise, tachyonic instabilities associated with a negative 
effective squared mass develop and prevent the existence 
of a static star configuration. This problem affects in 
particular the highly relativistic constant energy density 
stars. This is why we have used a polytropic equation 
of state, which is has the additional advantage to be a 
better approximation to a realistic neutron star. 

Our starting point is the action 
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mally coupled to the so-called Jordan metric g^u] R is 
the scalar curvature associated with g^^ and we have 
defined Mp^ = SttG. In contrast with Q and Q, we 
reexpress the model as a scalar-tensor theory in the so- 
called Einstein frame. The two formulations are of course 
equivalent (at least on the classical level), but the Ein- 
stein frame is useful to see that the behaviour of f{R) 
theories is quite analogous to chameleon-like scalar ten- 
sor theories. 

By introducing the scalar field = i/|il/pln/p 



and the metric 51^^ ~ Q, ^ g^^, , with 51 ^ = / ^ = 
^^/A/p] = exp[-2Q^/A/p], the action ^ be- 
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with the potential 
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which can be expressed in terms of by inverting the 
definition of as a function of R. 

The first and second derivatives of the potential V in 
terms of 6 are 
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The function /(i?) is severely constrained by obser- 
vations. Here we will consider the model suggested by 
Starobinsky 
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where Sm is the action for the matter, which is mini- 



For this model the potential ¥{(()) is shown in Fig. [T] 
Substituting this expression into ^ yields the minimum 
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of V corresponding to the asymptotic de Sitter solution 
_Roo = XocRo- It is convenient to express the parameter 
A in terms of x^^:. 



X = 
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where a prime denotes a derivative with respect to the 
radial coordinate, and one recognizes the total energy 
density in the brackets on the right hand side of the first 
equation. Instead of the angular component of Einstein's 
equations, one can use the conservation of the energy- 
momentum conservation, V aT't ~ 0, which reads 



In the strong gravity limit, corresponding to a; 3> 1, 
we have 
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where fJ^ = dfl/d(f). The Klein-Gordon equation for the 
scalar field reads. 



The curvature singularity thus corresponds to a finite 
value (j) = ioi the scalar field, as well as a finite value for 
the potential V = {X/2)MpRo. However, the derivative 
dV/ dcj) becomes infinite in the limit — > 0~ . 
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FIG. 1: Potential V (in units of MpRo) as a function of (in 
Planck units) for n = 1 and Xoo = 3.6. The lower black dot 
corresponds to the de-Sitter attractor while the upper-right 
dot shows the curvature singularity. 

For matter we consider a perfect fiuid characterized 
by the energy density p and pressure P, defined in the 
Jordan frame (i.e. with respect to the metric 5^1/). In 
the Einstein frame, the corresponding energy density and 
pressure are respectively p = fJ^p and P = fl'^P. Inside 
the star, (j) will be extremely close to zero and therefore, 
the quantities defined in the two frames will be numeri- 
cally indistinguishable. 

We consider a static and spherically symmetric metric 
(in the Einstein frame) 
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Introducing the radial function m(r) so that e^^ = 
1 — 2m/r, the time and radial components of Einstein's 
equations yield 
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In order to close the system of equations dH), ([7]), dS]) 
and ([9|) , one should finally specify the equation of state 
for the matter. In relativistic stars with constant en- 
ergy density were considered. The advantage is that the 
pressure can be easily determined analytically in general 
relativity (i.e. ignoring the backreaction of the scalar 
field) and is given by 
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where M is the mass of the star and r, its radius. How- 
ever, this implies that, in the innermost regions of the 
star, the quantity p~ 3P, i.e. the opposite of the trace of 
the energy-momentum tensor, becomes negative as soon 
as 
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In stars with a large region where p — 3P < 0, we were 
unable to find numerically solutions for the scalar field. 
This is not related to the presence of the curvature sin- 
gularity, since it occurs in chameleon type models as well 
p. The reason is the presence of instabilities, because 
the matter contribution to the effective squared mass is 
proportional to p — 3P and thus becomes ne gati ve. Note 
that a similar effect has been observed in for an- 
other class of scalar-tensor theories, where Infi depends 
quadratically on the scalar field . 

The equation of state deep inside a neutron star is not 
known but most realistic equations of state satisfy the 
condition p — 3P > 0. This is why we have chosen for the 
present analysis a polytropic equation of state, given by 



p{h) ^TTiB [n + K — , P{h) = KruB — , (10) 
no / no 



with ruB = 1.66 x 10"^^ kg, no = 0.1 fm~' and K = 0.1. 
In Fig. 2, the radial profile of the energy density and 
pressure, together with p — 3P, is plotted for a star with 
the central particle number density fic — 0.4 fm~^. This 
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FIG. 2: Energy density p (upper solid line), pressure P (lower 
solid line) and the combination p — 3P (dashed line), in units 
of the central density pc, as functions of the radial coordinate 
r (in units of Mpp^^^^). 
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corresponds to |$*| ~ 0.25, which is more than twice 
higher than the value reached in Q. 

We have integrated numerically the system of equa- 
tions (jUIll), together with (fTO|) . to obtain the profile of the 
scalar field inside and outside the star in the f{R) model 
(|4]) with 71 1. A key parameter is the ratio between the 
energy density at infinity (i.e. the cosmological energy 
density) and the energy density at the center of the star, 
which can be parametrized by vq = MpRo/ j5c- Realistic 
values of this parameter arc extremely small and are nu- 
merically challenging to explore, because the scalar field 
value at the center is proportional to Wq, as we will see 
below. The smallest value we considered is vq = 10~^, 
which is much higher than the realistic one. However 
we believe that the situation is qualitatively similar for 
smaller values. Instead of the shooting method used in 
Q , we have resorted to a relaxation method. 

As is clear from Fig. [31 the scalar field tends to inter- 
polate between an extremely high density regime, inside 
the star, and a very low density regime, outside the star. 
In the very high density regime, the scalar field is, nu- 
merically, very close to the singularity. This behaviour is 
quite analogous to that of the "chameleon" mechanism 
, as we will discuss in more details in Q . 

In order to understand intuitively our results, it is use- 
ful to introduce the effective potential 



(11) 



which is obtained by integrating, with respect to 0, the 
terms between brackets on the right hand side of the 
Klein-Gordon equation ^ (and we have set the constant 
of integration to zero) . The local minimum of the effec- 
tive potential, which we will denote (/'min(''), is charac- 
terized, if it exists, by 



d(j) d(h Mp' 



(p - 3P) = 0. 



FIG. 3: Profile of the scalar field (j> (in Planck units), shown 
by solid (blue) line, as a function of the radius (in units of 
Mppc^^'^), for the model ((U with n = 1, Xac ~ 3.6 and 
no ~ 10^*. The value (pmin for the minimum of the effective 
potential is plotted by dashed (gray) line. 



At high curvature {x ^ 1), dV/dcj) w AIpRq x/\/6, ac- 
cording to (O, and the minimum corresponds to 

P-3P 

which is inversely proportional to fo (the scalar field 
value (/)min IS thus proportional to Vq). This minimum 
exists only if the matter term p — 3P is positive. This 
implies that, for a very compact constant density star, 
the effective potential in the central layers of the star 
docs not have any mimimum. 

One can also define a corresponding effective squared 
mass as 
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which can become negative, as discussed earlier, if the 
second term dominates and p — 3P < 0. In the regime 
X ^ 1, if the minimum (|12p exists, one finds by evaluat- 
ing at the minimum, using ([5|), 
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In view of our results, which seem contradictory with 
the conclusions of and Q, it is worth reexamining 
in detail their arguments against the existence of very 
relativistic stars. Frolov's argument is based on the non- 
relativistic limit of the Klein-Gordon equation which can 
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be written in the form 



for t/i = /j=, — 1 w ~2Q4)/Mp. If one neglects the sec- 
ond term on the right hand side, ip then satisfies an 
equation analogous to the Poisson equation A$ = AnGp 
and is thus related to the gravitational potential $ by 
i' = "000 — 2$/3. Since tpao is small, this suggests that 
the curvature singularity can be easily reached with a 
sufficiently strong gravitational potential. However, in- 
side the star, the energy density is much higher than at 
infinity and this corresponds precisely to the high curva- 
ture regime, where dV/d(f> becomes very high and cannot 
be neglected. 




FIG. 4: Effective potential U{4>) = —V{(j)) in which the scalar 
field (j) "moves" under the action of the force in the picture 
of the classical mechanic analogy. 

One can also understand the scalar field profile, as in 
, in terms of a kinematic analogy with the motion of a 
particle in a potential U {(f) which is submitted to a force 
T = —{8ttG/3){p — 3P), where the radius r is consid- 
ered to be a "fictivc" time parameter (see Fig. The 
radial profile of the scalar field can then interpreted as 
the "motion" of the particle from r = to r = oo. At 
r = 0, the particle is at rest since 0' = 0, but is submitted 
to a leftward force that makes it roll uphill. The initial 
value for must be such that the scalar field reaches 
the top of the potential with vanishing velocity. We find 
a stronger gravitational potential in more compact stars 
with high energy density. The amplitude of then tends 
to become stronger, but there does not seem to be any 
critical value beyond which the particle would necessar- 
ily fall into the singularity or overshoot the top of the 
hill. As the density of the star increases, the minimum 
(f'min of the effective potential becomes closer to the sin- 
gularity 4> = 0, according to (fT2|) . One thus expects that 



the initial value of will also be closer and closer to the 
singularity = 0. This is confirmed by our numerical 
analysis, even if it becomes more and more difficult nu- 
merically to distinguish the extremely small value of (f> 
from zero. 

Moreover, as the numerical solution plotted in Fig. [3] 
shows explicitly, the "motion" of the particle can be sub- 
tle. In this case, it first starts to approach the singular- 
ity, it then stops and moves and the opposite direction 
towards the top of the potential. This complicated "mo- 
tion" is a direct consequence of the non-monotonous de- 
pendence of p — 3P, which reaches its maximum value at 
some non-zero radius of the star, as can be seen in Fig. [5] 

To conclude, our work shows that it is possible to con- 
struct highly relativistic stars in /(i?) theories, despite 
recent indications to the contrary. Numerically, the task 
can be challenging as the scalar field value is extremely 
close to the singularity in the center of the star. We have 
also studied relativistic stars in chameleon-like models 
and found a very similar behavior Q. In both cases, 
a crucial requirement for obtaining a static configura- 
tion is that the equation of state satisfies the condition 
p — BP > in most of the star. Otherwise, the effec- 
tive squared mass of the scalar field becomes negative 
and the associated instabilities prevent the construction 
of a static configuration. This is problematic for con- 
stant energy density stars, which have been used in sev- 
eral recent works because of their analytical simplicity, 
since p — 3P becomes negative in the central part of very 
massive stars of this type. This might explain why, be- 
yond some critical value for the gravitational field, no 
configuration was found in for chamelon-like models. 
Realistic neutron stars do not, however, suffer from this 
problem. Moreover, once the scalar field configuration is 
stable, the (hydrodynamical) stability of the stars should 
be very similar to that of their general relativistic coun- 
terparts because the backreaction of the scalar field on 
the star is extremely small. An open question, which is 
far beyond the scope of this work, is whether the static 
configurations described here can be reached dynamically 
during the collapse of a very massive star into a neutron 
star. It would also be interesting to study whether and 
how black holes would form during the collapse of mat- 
ter, i.e. such has been done for instance in [l^ in the 
context of scalar-tensor theories without potential. 
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